Interference Coordination: Random Clustering and 

Adaptive Limited Feedback 

Salam Akoum and Robert W. Heath, Jr. 



O 

(N 

> 

O 

O: 

m 
(N 



> 

o 
d 

(N 



X: 



Abstract — Interference coordination improves data rates and 
reduces outages in cellular networks. Accurately evaluating the 
gains of coordination, however, is contingent upon using a 
network topology that models realistic cellular deployments. In 
this paper, we model the base stations locations as a Poisson 
point process to provide a better analytical assessment of the 
performance of coordination. Since interference coordination 
is only feasible within clusters of limited size, we consider a 
random clustering process where cluster stations are located 
according to a random point process and groups of base stations 
associated with the same cluster coordinate. We assume channel 
knowledge is exchanged among coordinating base stations, and 
we analyze the performance of interference coordination when 
channel knowledge at the transmitters is either perfect or ac- 
quired through limited feedback. We apply intercell interference 
nulling (ICIN) to coordinate interference inside the clusters. The 
feasibility of ICIN depends on the number of antennas at the 
base stations. Using tools from stochastic geometry, we derive the 
probability of coverage and the average rate for a typical mobile 
user. We show that the average cluster size can be optimized as 
a function of the number of antennas to maximize the gains 
of ICIN. To minimize the mean loss in rate due to limited 
feedback, we propose an adaptive feedback allocation strategy 
at the mobile users. We show that adapting the bit allocation 
as a function of the signals' strength increases the achievable 
rate with limited feedback, compared to equal bit partitioning. 
Finally, we illustrate how this analysis can help solve network 
design problems such as identifying regions where coordination 
provides gains based on average cluster size, number of antennas, 
and number of feedback bits. 



I. Introduction 

Coordination can mitigate interference and increase data 
rates in cellular systems [jj. Complete coordination between 
all the base stations in the network, however, is not feasible 
121, lis). For inter-base station overhead to be affordable, 
only groups of base stations coordinate, forming coordination 
clusters ||2l, ||4l, Q. To quantify the gains from coordination, 
an accurate model of the network topology and the relative 
locations of the mobile users and the base stations needs to 
be considered H . Much prior analytical work on interference 
coordination used oversimplified network models and reported 
coordination gains that did not materialize in practical cellular 
deployment scenarios ||6l. This paper addresses this issue 
by considering a point process model for deployment and 
clustering of base stations. 

Most of the literature on interference coordination H], lH, 
IEIj Q-fSl considered fixed cellular network architectures 
such as the Wyner model or the hexagonal grid. The Wyner 

S. Akoum and R. Heath are with The University of Texas at Austin, 
2501 speedway stop C0806, Austin, TX 78712 USA (e-mail: salam.akoum, 
rheath@mail.utexas.edu). This work was funded in part by Huawei Technolo- 
gies. 



model is used to derive information theoretic bounds on the 
performance of multicell cooperation HI, ID, [Til; it does 
not account for the random locations of the users inside the 
cells. Hexagonal grid models, although reasonably successful 
in studying cellular networks, generally rely on extensive 
Monte Carlo simulations to gain insight into effective system 
design parameters 121, ifTSl . Using random models for the base 
stations locations yields, under fairly simple assumptions, ana- 
lytical characterizations of outage and capacity, and provides a 
good approximation for the performance of actual base stations 
deployment |[T6l . 

We consider randomly deployed base stations. To form coor- 
dinating base station clusters, we propose a random clustering 
model, in which cluster stations are randomly deployed in 
the plane, and base stations connect to their geographically 
closest cluster station. The random clustering model builds on 
the analytical appeal of the random network deployment. It 
mirrors the connection of base transceiver stations to their 
base station controllers in current cellular systems ITtI . A 
similar hierarchical structure was used in fTsl to model traffic 
and requests for communication in wired telecommunication 
networks. A regular lattice clustering model of randomly 
deployed base stations was considered in |T9l to derive the 
asymptotic outage performance of interference coordination 
as a function of the location of the user inside the fixed grid 
and the scattering model. 

To achieve coordination gains, coordinating base stations 
exchange channel state information (CSI) on the backhaul. 
Each base station designs its beamforming vector to transmit 
exclusively to the users in its own cell, while exchanging CSI 
of the served users with the cooperating base stations. Ex- 
amples of interference coordination strategies in the literature 
El-IIOl, im, ini include iterative strategies such as MMSE 
estimation beamforming flOl and non-iterative strategies such 
as intercell interference nulling (ICIN) 13, El, d. ICIN is 
a cooperative beamforming strategy where each base station 
transmits in the null space of its interference channels. 

The gains from coordination depend on the quality of CSI 
available to the base stations. In frequency division duplex 
systems, the CSI of the estimated channels at the mobile users 
can be made available at the base stations through feedback. 
The feedback channel possesses finite bandwidth and thus 
limited feedback techniques 1201 are employed. Most of the 
literature on multicell coordination Q, ITTl . lT4l assumes full 
CSI at the transmitters, neglecting the performance loss due 
to quantization. For ICIN with limited feedback, each mobile 
user quantizes and feeds back the CSI of the desired and 
interfering channels. This is in contrast with single-cell limited 
feedback , where only the CSI of the desired channel is fed 



back Il2n . Dividing the feedback budget among the desired 
and interfering channels has been proposed as a reasonable 
approach for ICIN with limited feedback IHl, lH. The 
allocation of bits to channels at the receiver can be done 
equally or adaptively, taking into account the location of the 
user and the strength of the channel being quantized f 131 . 

In this paper, we analyze the performance gains of inter- 
ference coordination in a clustered network model where the 
base station locations are distributed as an independent homo- 
geneous Poisson point process (PPP) and the clustering of the 
base stations is done through overlaying another independent 
homogeneous PPP. Base stations belonging to the same cluster 
coordinate via ICIN. We account for the feasibility of ICIN 
given the number of interferers in each cluster and the number 
of antennas at each base station. We distinguish between 
two cases, first assuming that ICIN is always feasible, and 
later considering a preset number of antennas at each base 
station and only applying ICIN when feasible. Single-cell 
beamforming is applied otherwise. Using tools from stochastic 
geometry, we derive a bound on the coverage and average rate 
performance of the proposed clustered coordination model for 
both cases. For the case of limited feedback CSI, we derive 
a bound on the mean loss in rate. We assume random vector 
quantization (RVQ) at the mobile users for analytical reasons. 
To minimize the mean loss in rate due to quantization, we 
consider an adaptive bit allocation algorithm and we derive 
closed form expressions to determine the number of bits to 
allocate to each channel based on its strength. 

The contributions of the paper are summarized as follows. 

• We propose a random hierarchical clustering model based 
on point process deployment of cluster stations and 
base stations. We assume that the base stations connect 
to their geographically closed cluster stations to form 
coordination clusters. 

• We derive bounds on the coverage and average rate 
performance of ICIN, for the proposed clustered model, 
as a function of the number of interferers inside each 
cluster, the number of antennas at each base station, and 
the average cluster size. 

• We analyze the feasibility of ICIN taking into account 
the random number of interferers inside each cluster and 
the number of antennas at the base stations. We propose 
a thresholding policy where ICIN is applied only when 
feasible, and single-cell beamforming is applied other- 
wise. We derive bounds on the probability of coverage 
and average rate performance of the thresholding policy. 
We show that there is an optimal cluster size to achieve 
maximum coordination gains. 

• We derive an upper bound on the mean loss in rate due to 
limited feedback using ICIN. The bound is a function of 
the number of feedback bits allocated to each channel, the 
number of interferers, and the relative channel strengths 
at the mobile user We derive closed form expressions for 
allocating the feedback bits to the desired and interfering 
channels, as a function of their relative strengths at the 
mobile user, and taking into account the relative strength 
of the inter-cluster and intra-cluster interference. 
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Fig. 1. The hierarchical cellular model considered. The Voronoi tessellation 
of the plane formed by the cells of lib is shown in blue. The Voronoi 
tessellation formed by EIc is shown in red. The typical mobile user is shown 
in a green circle. The cooperating base stations for the cluster of interest are 
shown in rectangles. 



Throughout the paper, we use the following notation, a is 
a column vector, A is a matrix, a is a scalar. A* denotes 
the conjugate transpose and A~^ the matrix inverse. The left 
pseudo-inverse of A with linearly independent columns is 
= (A*A)~^ A*. E{.} denotes expectation. |.| denotes the 
cardinality of a set. 

II. System Description and Assumptions 

Consider the network model shown in Figure [T] The base 
stations are represented by an independent homogeneous PPP 
lib of density Ab. The mobile users are located according 
to an independent point process IIu. The cells of the base 
stations form a Voronoi tessellation of the plane with respect 
to the process lib, ™d the users are connected to their closest 
base station. We assume that the density of the mobile process 
is sufficiently large such that all the base stations are active. 
Each base station is equipped with A^t antennas and serves 
one single antenna mobile user in each cell using intra-cell 
time division multiple access (TDMA). While TDMA is not 
necessarily the best option for a multiple-input-single-output 
(MISO) transmission strategy, it is a common assumption 
made in the multicell cooperation literature due to the tractabil- 
ity of single-user transmission HI, lfT3l . J?], |fT9l . 

To model clustering, we overlay another independent homo- 
geneous PPP He on the 2-D plane with intensity Ac < Ab. We 
denote the cluster centers of this PPP by Ck- The cluster cells 
form a Voronoi tessellation of the plane with respect to the 
process He. The base stations of the lib process located in the 
same Voronoi cell l/t'^fc) of a cluster center Ck coordinate via 
ICIN. The association of the base stations to the cluster base 
stations follows similar to the association of the mobile users 
to their closest base station in the basic cellular model. The 
cluster base stations can be seen as central processors to which 
base stations in the same cluster are connected via backhaul. 
The backhaul is assumed error and delay free. 

Based on the stationarity of the Poisson process 1231 . we 
consider the performance for a typical user uq served by 
a base station bo inside a cluster ^^'^''^(bo). The channel 
corresponding to the desired signal between bp and uo is 
denoted by hp G C^^^'. The interfering channels from the 
^-th base station to uq are denoted by ^ e C'^''^\ The 
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desired and interfering channels are modeled according to the 
Rayleigh fading model, where each entry is independently and 
identically distributed as a unit variance zero mean complex 
Gaussian. The symbol transmitted from the fc-th base station 
is denoted by .s^, where E[|sfcp] = Es and there is no 
power control. While power control is an important topic for 
practical networks, an exhaustive investigation of the topic 
would require a detailed analysis of a distributed power control 
algorithm among the base stations. Another way of imple- 
menting power control independently at each base station is 
to perform channel inversion and assume fixed received power 
ll24l . Such a strategy leads to better outage performance for 
full CSI |[T9l . In the absence of interference coordination, each 
user is subject to interference from all other base stations in 
lib, each transmitting with energy Eg. The path-loss incurred 
by the desired signal is given by L(ro) = (1 + fp)" where 
ro is the distance between uq and bo, and a is the path-loss 
exponent. For the interfering signals from the ^-th base station 
to uo, the path-loss is given by L{roj) = {l+roj)" ,roj > ro- 
Most of the literature on analysis of random spatial models 
|[T6l . |fT9l consider the L{r) = r" path-loss model; this 
model, however, has a singularity at zero and is inaccurate 
for small distances. The received signal powers of the desired 
and interfering signals are then given by 70 = Es/L(ro) and 
7o.f = Es/L(ro.£). Using a narrowband flat-fading model, the 
baseband discrete-time input-output relation for uq is given by 



2/0 = y/^Kioso 



/Tu7g5,/«Sf + f'o, (1) 



where yo is the received signal at uq, the vector fo G C^'^^ 
is the unit-norm beamforming vector at bo, and G C^'^^ 
is the unit-norm beamforming vector at the £-th base station. 
The scalar vo denotes the additive white Gaussian noise at the 
receiver with variance a^. 



The signal-to-interference-plus-noise ratio (SINR) at uq is 
given by 



SINRo 



TolhSfoP 



+ Eb,enb/boTo,«|g5/« 



(2) 



For a given SINR threshold T, a performance metric of interest 
is the probability of coverage. 



p,(Ab,Ac,a,r) ==P[SINRo >T] 



(3) 



The probability of coverage is the probability that a randomly 
chosen user in the 2-D plane has a target SINR greater than T. 
Such a condition is required in practice when a given constant 
bit-rate corresponding to a particular coding scheme needs to 
be sustained. 



Another performance metric of interest is the average rate, 

T(Ab,Ac,a) =E{log2(l + SINRo)}. (4) 

To express the average rate, we assume that adaptive modula- 
tion/coding is used so that the users achieve a Shannon bound 
for their rate. 



III. Per Cluster Interference Coordination 

We propose interference coordination within each cluster 
based on zero forcing (ZF) intercell interference nulling. We 
assume there are N interfering base stations in the cluster of 
interest. Since is a random variable that depends on Ac 
and Ab, the number of antennas at the base stations needs to 
be sufficiently high (A^t > N) for ICIN to be feasible inside 
each cluster. We assume for the first part of this paper that 
A't grows with N such that A't — iV = d^Vt, where dAr^ are 
the extra degrees of freedom at each base station, used to 
beamform to the desired user ||251 . 1261 . This assumption is 
relaxed to a preset constant A't in Section|VT] Each base station 
in the cluster has knowledge not only of the channel to its 
intended receiver, but also of the interference channels towards 
the receivers in the same cluster Without loss of generality, for 
the case of the typical user, uq estimates its desired channel 
ho and the interfering channels go.^, £ = !,■■■ ,N, and feeds 
the information back to its base station bo. The base stations 
in the coordination cluster then exchange the interfering CSI 
on the backhaul, so that bo has knowledge of ho and the 
interference channels towards other receivers in the cluster 
gf,oJ=l,--- ,N, 

We assume that perfect CSI is available at the receivers. 
Accounting for imperfections due to channels downlink train- 
ing further degrades the achievable rate with imperfect CSI, 
and is a subject of future investigation. We consider first the 
ideal case where full CSI is available at the transmitters, we 
then discuss finite rate feedback. 

A. ICIN with Perfect CSI 

With perfect CSI at the transmitters, each base station 
designs its beamforming vector such that it is the normalized 
projection of the desired channel direction onto the null space 
of the matrix of interfering channel directions. For bo, fo is the 
normalized projection of ho = ho/||ho|| onto the null space 
of G/,0 = [gi,o • • ■ gefi ■ ■ ■ gN,o], with gift = Se,o/\\ge,o\\, 



fo = 







|ho 






|ho 



(5) 



The SINRo with per cluster ICIN, is given by 

7o|hSfor 



SINRic = 



Eb,enb/v(-o)(bo) 7o,£|g(*,/«P + cr^ 
(l + ro)-"|hSfop 



(6) 



/out -h 1/SNR ' 

where /out is the inter-cluster interference caused by all the 
base stations outside the cluster of interest V^'^"^ho), denoted 
by (ho). The signal-to-noise-ratio (SNR) here is de- 

fined in terms of the transmit signal power, 1/SNR = (T^/Eg. 

B. ICIN with Limited Feedback CSI 

With limited feedback, the estimated channel directions ho 
and go,^ are quantized to the unit-norm vectors given by 
ho and goj, respectively, at uo. These quantized channel 
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directions are then fed back to the base stations using a fixed 
feedback budget of Btot bits. The base stations use these quan- 
tized channel directions to design the beamforming vector We 
assume that the base stations have perfect knowledge of the 
SNR at the receivers, independently of the channel directions. 
Perfect non-quantized knowledge of the received SNR at the 
base stations is a common assumption in the literature on 
multi-user and multicell limited feedback ifTSl . 1271 . 

We consider separate quantization at the receiver of the 
desired and interfering signals. Joint quantization was con- 
sidered in ||28l ; it requires a large storage space at the mobile 
users. Each user divides Btot among the desired and interfering 
channels such that Bo and Bo,f are used to quantize hg 
and go.f respectively, and Bp + X^fci ^oj ~ Btot bits. We 
ignore delays on the feedback channel as well as delays on 
the backhaul link between the base stations belonging to the 
same cluster The beamforming vector fo is the normalized 
projection of the quantized channel ho onto the nuU space of 
the matrix G/,o = [gi,o ■ • • gi^o ■ ■ ■ gN,o], 



fo 





~ G/,oG|_(,^ 


1 ho 




- G/,oGJ_(,^ 


1 ho 



(7) 



SINRo with limited feedback and per cluster ICIN, is given 

by 



SINRic 



(l + ro)-"|hSfop 



where /^cs = Eb,Gy('^o)(bo)/bo (1 + ''o,^) "\So/e? is the 
residual intra-cluster interference due to quantization. Interfer- 
ence signals are not nulled out, since gjj ^f^ ^ but gp ^f^ = 0. 

Tout = T,hcenjv(-o)(ho\i^ + '^o,ty"\go/e\^ is the inter- 
cluster interference assuming limited feedback beamforming. 

We derive the average rate with perfect CSI in Section |IV] 
We then bound the mean loss in rate with limited feedback 
for equal bit and adaptive bit partitioning among the desired 
and interfering channels in Section [V] 

IV. Performance Evaluation with Perfect CSI 

In this section, we analyze the outage performance of 
per cluster ICIN. We then use the probability of coverage 
expression to derive a bound on the average rate. 

A. Probability of Coverage Analysis 

The probability of coverage is the complementary cumula- 
tive distribution (CCDF) of SINRic, 



The aggregate inter-cluster interference power /out is a 
function of the cluster size. While a numerical fit of the 
distribution of the cluster size is possible 1291 , an expression 
for the exact size of the Voronoi region 1/(^=0) (bo) is hard to 
compute. Hence, we bound the area of the cluster of interest 
by the area of the maximal disk inscribed in the cluster, 
and centered at Co. Let 7'i„ designate the radius of Bi„. r„^ is 
Rayleigh distributed with CCDF P[r,„ > r] = cxp(-47rAcr2), 
Il30i . The interference power is then upper bounded by 



(9) 



/out = ^ ^o,e\So,Af < X] lo,i\So,A 
nb/i'(^)(bo) nb/B„ 

< ^ 'yo,e\so,efe\^, 
nb/B,„{uo) 

where Bm(uo) denotes the disk centered at the typical mobile 
user uo of radius Vi^ — ro — ri, where ri is the distance from bo 
to Co, Rayleigh distributed with probability density function 
(PDF) fri{r) ~ 27rAc7' cxp(— TrAcr^). The interference field 
power is computed at uo located at a distance < r < 
from the center of the inscribed disk Bm corresponding to 
the cluster-cell center The exclusion distance to the nearest 
interferer from the mobile user is asymmetric since the dis- 
tance from the mobile user to the closest edge of the disk 
is smaller than its distance from the furthest edge. To avoid 
the dependence on the location of the interference field, we 
pursue an upper bound on the interference power We consider 
the interference contribution in a disk of radius 7'i„ — ri — vq 
centered at uq and denoted by Bni(uo). The exclusion areas 
are such that Bm(uo) C B^ C y(^°)(bo). The upper bound 
on the aggregate inter-cluster interference power results in a 
lower bound on the probability of coverage. 

Result 1: The probability of coverage with per cluster in- 
terference coordination is lower bounded by 

Pc.ic(Ab,Ac,a,T) > W fr{ro) /r„(rm)x 
r e-^-^^'Cj (2i.L(r-o)T.)M:^^!^dsdr.odr„ 



(10) 



where the Laplace transform of the desired signal Ibpfop 

r[djv„i], is 



Ch{s) 



Nt-N. 



(s + l)°"t 

The PDF of tq, fr{r) = 27rAb?'exp(— TrAbr^), follows from 
the null probability of a 2-D PPP with respect to lib. The 
PDF of the radius r,„ is fr^{r) = SttAcT cxp(— 47rAcr^). The 
Laplace transform of the interference Cj^ is 



Pc(Ab, Ac, Q,T) = P[SINRic>T] 

= P[|hSfo|' >r/(r-o)(/out + l/SNR) 



(8) 



The desired effective channel power |hofoP is Gamma- 
distributed with parameters djy^ and 1, lligfop ~ r[d7Vj,l]. 
This follows from the projection of the independent normal 
Gaussian isotropic vector ho onto the null space of the 
interfering channels of dimension dAr^ = Nt — N ET\ . Il25l . 



exp — 27r 



(11) 



Ab 

(l + r)^-" 

Q - 1 



(l + r)2 



S2F1 



2 

1, 1 



s 2F1 
1 



-.2 - 



a 



1, 1 - 
1 

a ' 



r(c 



where 2Fi(a,6,c,z) - r(b)r(c-b) Jo 
the Gauss hypergeometric function. 



2 2 

- -,2 ,- 

a a 

"1 t''"^(l-t) 



-(1 + r)- 



(l-tz)" 



denotes 
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Proof: See Appendix lAl ■ 
The tightness of the upper bound on the probability of outage 
is investigated in Section IVIII The probability of coverage 
depends on the density of interfering base stations, the density 
of the clusters, the path-loss exponent and the number of extra 
spatial dimensions available at the transmitters d^Vf It is an 
increasing function of the ratio of the density of interfering 
base stations to the density of the cluster base stations Ab/Ac, 
which is also the average number of base stations per cluster 
It is also an increasing function of cIat^. As A^t increases, djVj 
increases, and the signal power at the mobile users increases, 
which incurs an increase in the SINR, in addition to the 
decrease in the interference term due to ICIN per cluster 

B. Average Rate 

The average rate follows from the probabiUty of coverage 
analysis as follows 



E{log2 (1 + SINRo)} 



f 
Jo 



' [SINRo > 



log(2) 



'-dt. (12) 



Inserting (fTOl i into ( fT2] i gives a lower bound on the average 
rate of the typical user. 

Result 2: The throughput of the typical user uq, averaged 
over the spatial realizations of the point processes lib and He, 
and the channel fading distributions is bounded by 



f 

Jo 

r 

Jo 



¥ [SINRic > - l] 

Pc.ic (Ab, Ac,a,t)) 

(l + i;)log(2) 
plb.^ {Ab,Ac,a,i.) 



dt 



dv 



dv, 



(13) 



(l + i))log(2) 

where p^^-^^ (Ab, Ac, a, v) is the lower bound on the probability 
of coverage pc.ic (Ab, Ac, a, v) on the right hand side of ( fTOl i. 

The computation of Tic requires an additional numerical 
integration over Pc,ic- As in the case of coverage. Tic depends 
on the average cluster size. It increases with increasing Jat^ 
and increasing Ab/Ac, as illustrated in Section [VTIl 



V. Impact of Limited Feedback on Average Rate 

The performance of ICIN depends on the availability of 
accurate channel state information at the cooperating base 
stations. With limited feedback, the quantized CSI is fed back 
from the receiver to the transmitter. The feedback resources 
at each mobile user are partitioned among the desired and 
interfering channels. The bit allocation strategy between the 
quantized channels is hence important to minimize the mean 
loss in rate due to limited feedback. 

We compute an upper bound on the mean loss in rate 
using limited feedback. We use separate codebooks to quantize 
each desired and interfering channel. We consider random 
vector quantization in which each of the quantization vectors 
is independently chosen from the isotropic distribution on the 
A^t dimensional unit sphere ISTl . RVQ is chosen because it 
is amenable to analysis and yields tractable bounds on the 
mean loss in rate. Furthermore, optimum codebooks are not 
yet known for multicell cooperative transmission llT3l . 



We define the mean loss in rate due to limited feedback as 



= E<^ log2 (1 + SINRic) \ - E<^ log2 (l + SINR; 



(14) 

where Tic is the average rate with limited feedbaclfl 

The SINRs with perfect and limited feedback CSI are such 
that SINRic > sTnRic. Subsequently log2(l + SINRic) - log2(l + 
SINRic) < log2(SINRic) - log2(Sll\IRic) and 

Afic < E{log2(SINRic) - log2(SlNRic)}. (15) 
We denote the upper bound on ATic by At^"*^ given by 

lowt + cSp + -frcs 



Af.".''=E^ log2(SM 



|h2fo 



+ M log2 



/out + sRr 



where AtcIcs denotes the mean loss from quantizing the 
desired channel and ATint is the mean loss from quantizing 
the interference channels. 

To derive the contribution of the desired signal quantization, 
ATdos, we use a lower bound on the desired signal power with 
limited feedback, lUS 

+ E [log2 (||ho|p|h5fop)] . (16) 
An upper bound on ATdcs then follows 



Ardos < log2(e) Ejv 



Aft - 1 t-r V A^t - 1 



(17) 



since E|log2 (llholPlhSfop) I = e| loga (||ho|p|hSfop 

The Beta function is defined as /3(a, b) = - tf-^dt. 

The effective interference power from each base station 
outside the cluster is a function of the interference channels 
towards the typical user, and the quantization vector f^. The 
quantization vector is independent of gg ^ and is computed 
using RVQ. The interference power from each base station 
is thus distributed as |go/fP ~ r[l,l],€ G I,-- - ,N, 
131:1. We rewrite ATi^t as 



Arint = E<j log2 ( /out + + has 



1 

SNR, 

and we focus hereafter on the effect of the residual intra-cluster 
interference I-^^s- 

The statistics of /res and the bound on the mean loss in 
rate depends on the number of interfering base stations inside 
the cluster N, and the strategy of allocating bits among the 
desired and interfering channels. We first consider an equal bit 
allocation strategy, independent of the channels' strength, we 

'e/ log2 ^1 + SINRic^ J- is an upper bound on the achievable rate with 

limited feedback |32| . This bound is approached when the mobile user has 
perfect knowledge of the coupHng coefficients between the beamforming 
vectors used at the base stations and the channel vectors 1321 . 
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then optimize the bit allocation to minimize the mean loss in 
rate. 



A. Equal Bit Allocation 

One option for bit allocation is to divide Btot almost equally 
between the interfering channels and the desired channel, that 
is Bo,£ = LBtot/(A^ + 1)J and then set Bq = Bft - NBo,e- 
This is a suboptimal strategy because it does not account for 
the difference in the path-loss between the various interferers, 
and hence the difference in the importance of the interfering 
signals. It gives, however, a slight bias for quantizing the 
desired channel, when the number of bits is not a multiple 
of the number of base stations in the cluster. In this section, 
we compute a closed form expression for the bound on the 
mean loss in rate, and we use it as a stepping stone to later 
optimize the bit allocation. 

The bound on the mean loss in rate is rewritten as 



< log2(e) Ejv 



Nt - 1 



E<^ log2 (/out + 1/SNR) 



; (1 + ro,e)-°'\s*Q,M + lout + 1/SNR 



As At^"'' is averaged over all the spatial realizations in the 
network, with different number of interferers N, we first derive 
the probability mass function (PMF) of in a typical cluster 
The PMF of TV depends on the cluster size. It is the distribution 
of the number of points inside a cluster area, given that a 
randomly chosen point is located in that area. 

Lemma 1: The probability mass function of the number of 
interferers inside the cluster of interest is given by 



P[iV 



3.5^-^r(n + 4.5)(Ab/Ac)" 
r(4.5)n!(Ab/Ac + 3.5)"+4 



(18) 



Proof: See Appendix IbI 



For the contribution of the inter-cluster interference /out and 
consequently in AtJ"'', we approximate the distribution 
of I,.^ with the Gamma distribution, using second-order mo- 
ment matching Ii33ll . The Gamma approximation yields more 
tractable expressions for E{log(/r„)} than the challenging to 
compute density and Laplace characterizations of 

Remark 1: (Gamma Distribution Second Order Moment 
Matching) For the second-order moment matching of the 
Gamma random variable, consider a random variable x with 
finite first E{.t} and second order E{a;^} moments, and 
variance var(a;) = E{a;^} — (E{.t})^. The Gamma distribution 
r[fc, with the same first and second order moments as x has 



parameters k 



, and 



var(a:) 



For the residual intra-cluster interference, we invoke 
Jensen's inequality on the first term on the right hand side 
of ( fTSl l to derive a closed form expression for E{/rcs} 



(f) , 



1=1 

N 



where (a) follows from using RVQ for quantization 
^{\SoJe\'} = ^2SM/?(2«o,.^_ii_)^ ini, |3i|. (b) fol- 
lows from Stirling's approximation on the Beta function. 

Result 3: The mean loss in rate from limited feedback with 
ICIN, with equal bit allocation, is bounded by 



A?ic <log2(e)^r 



A^t 



A^t - 1 



2Nt ~ 1\ - 1 °tot I 1 

-J^^P ^'•+i^"*-iPiv(n)E{(l + ro,i)-"} 



(1 + r) 



^'^ fro (T-o)drdro, 



- log2(e)t/. (fc) - log2{e) + log2 ( 1/SNR + ke 

OO 

where 

and /out ~ r[fc,(?] via second-order moment matching 1331 . 
with parameters k and such that Ej/j.^^^} = kO and 
var(/r„J = kO'^, with 

E{/,„,} = 2^Ab E,„,,„ -^dr} and var(4„) ^ 

Proof: See Appendix ICl ■ 
The mean loss in rate Aric decreases as the total number of 
bits Btot increases. Atjc increases, however, with the average 
cluster size Ab/Ac, given a fixed feedback budget Btot- A 
higher Ab / Ac implies more base stations per cluster and hence 
more interfering channels to quantize. As the number of 
antennas at the base stations increases with the number of 
base stations per cluster for ICIN to be feasible, the number 
of bits allocated per antenna for each quantized channel also 
decreases. 

B. Adaptive Bit Allocation 

As the equal bit allocation limited feedback strategy results 
in a considerable decrease in the ICIN performance, as will 
be illustrated in Section [VTIl in this section we optimize the 
bit allocation to minimize the mean loss in achievable rate of 
ICIN. We adapt the number of bits allocated to the desired and 
interfering channels as a function of their signal strength at the 
typical user Since each spatial realization results in a different 
number of intra-cluster interferers N, the optimization is done 
per spatial realization, as a function of TV. 

We denote the total number of bits allocated to quantizing 
the interfering channels by B^ = Btot — Bq. Given B;, we 
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first derive Bo/,€ — I,-- - such that the contribution 
of the interfering channels towards the mean loss in rate is 
minimized. In other words, we aim at finding Bo,^ such that 
/res is minimized. The optimization problem is expressed as 



mm 



5^(i + ro,f)-"r 



2Nt 



Nt - 1 

s.t. EfllBo.f = and Bqj > 



2 «t-i 



(19) 



The bit allocations Bo,£,^ G I-- - ,iV are integer-valued. 
We solve a relaxation of the optimization problem assuming 
Bq/ are real valued. The solution is given by the arithmetic- 
geometric mean inequality. It is derived in |13| and restated 
here for completeness. 

Lemma 2: ( lfT3l Theorem 4]) The optimum number of bits 
assigned to the ^-th interferer, Bq ^, that minimizes ( fT9] l is 
given by 



for £ G /C and Bo,£ for i 
of interferers that satisfies 



(20) 



n«x(i+'-(M)-/i'ci 

: JC, where JC is the largest set 



log2 



(l+ro,^)-" 



< 



B, 



|/C|(iVt-l)- 



(21) 



JC denotes the set of effective interferers, such that \JC\ < N. 
Using the expression for Bqj from ( l20l i. the bound on the 
mean loss in rate per spatial realization, i.e. given N and ro,£ 
is expressed as a function of a single variable Bq = Btot — B;, 



Afi,(Bo) 



Nt 



Nt - 1 



2 "t-i 



■E<{ loga (/out + l/SNR) 

/2Aft-l\ Btot -Bp ,. . 

- log2 I r ( ) |/C|2 n (1 + ro,^)-"/l'=l 



+E{/out} + l/SNR). 



(22) 



To minimize the mean loss in rate, we minimize Aric(Bo) 
as a function of Bq such that Bq < Btot- We distinguish 
between two cases of interest, depending on the ratio of 
residual interference /res to inter-cluster interference /out. 

Dominant Inter-Cluster Interference: When E{/rcs} < 
E{/out} + 1/SNR, we use the low-SNR approximation log(l + 
x) w X. The bound on the mean loss in rate is written as 



Afic(Bo) S:ilog2(e)r 



Nt 



Nt-1 



2 "t-i 



- e| log2 (/out + 1/SNR) I + log2(E{/out} + 1/SNR) 
+ log2(e)^ 



Nt-l 



E{7out} + 1/SNR 



\IC\2 l'=l(«t-i) Yl (1 +r-o,«)"°'''l'^l 



ieic 



(23) 



The optimization is carried out for the terms in Bq. 

Result 4: The optimal value of Bq for E{/ros} < E{Jout} + 
1/SNR, given N and ro^ej^ 1 • • • iV, is 



Btot 



Nt\JC\ 



\K\ 



\K\ + 1 \ A^t 

+ log2(E{/out} + 1/SNR) 



Yn(i+^o,*)" 

1=1 
{Nt - 1)|/C| 
|/C| + 1 



{Nt-1)\K\ 



1/C|- 



(24) 



Proof: The proof is provided in Appendix [D] ■ 
Dominant Residual Intra-Cluster Interference When 

E{/res} > E{Jout} + 1/SNR, we use the high SNR approxi- 
mation log(l + x) Ri log(x) and we write ATic(Bo) as 



A?ic(Bo) «log2(e)r 



A^t 



Nt - 1 



-Iog2 



/ 2A^t -1 
V A^t - 1 



2 «t-i _ E<{ log2 (/out + 1/SNR) 
Bo - Btot 



1^1 n (1+^0/)"° 



/IKI 



\K.\{Nt-l) 
(25) 

Since the terms which are a function of Bq in dZSl l are convex 
with respect to Bq, the optimal value of Bq is found by taking 
the derivative of these terms with respect to Bp. 

Result 5: The optimal value of Bq for E{/res} > E{Jout} + 
1/SNR, given N, is 



Bo = (Art-l)log2 |/C|log2{e)r 



Aft 



A^t - 1 



(26) 



The optimization problems for both dominant inter-cluster and 
dominant intra-cluster interference are solved assuming non- 
integer bit values. Since the objective functions in (1231 ) and 
dZSl l are convex in Bq, we get the integer value of Bq by 
taking the floor or the ceiling of the value in ([M]) and ( l26b . 
respectively. 

An upper bound on mean loss in rate with adaptive bit 
allocations, averaged over all spatial realizations, is subse- 
quently given by replacing Bo by its optimal value in ( |23] ). 
distinguishing between the two dominant interference cases. 



VI. Performance Evaluation with Perfect CSI 
WITH Fixed iVt 

In the analysis so far, we assumed that iVt is a function of 
the number of interferers in the cluster and hence is a random 
variable. This assumption makes the analysis with perfect and 
imperfect CSI more tractable, and allows system designers to 
gauge how many antennas they should provide at the base 
stations for interference coordination to be most beneficial. 

In this section, since the number of antennas may be 
predetermined, we relax this assumption and we assume A't 
is fixed at each base station. As iV is a random variable that 
changes per spatial realization and can become greater than 
A^t, interference coordination may no longer be feasible per 
cluster. To overcome this issue, we propose a thresholding 
policy wherein if A^t < N, each base station beamforms to 
its own user and no interference coordination is performed. 
When Nt > N, ICIN is applied as in Section |IV] 

Let A be the event that the number of interferers in the 
cluster is less than Nt, then for a fixed Nt the probability of 
coverage is expressed as 
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' [SIRo > T\A]V [A] + P [SIRo > T\A'']F [A"] 
•[SIRo >TnA]+P[SIRo >rnA'=]. (27) 



The probability of coverage with per cluster coordination and 
Nt > N follows from the probability of coverage expression 
derived in Section |IV] with one difference being that the 
averaging is now done over N, given that A^t is fixed and the 
Laplace transform of the desired signal depends on iVt — N, 
now a random variable. Using the PMF of N in Section |V] 
the lower bound on the probability of coverage Pc,ic,Nt is 



Nt-l 

Pc,ic,iVt 



poo roo 

PN{n) / /r(r-o) / /im(rm)x 

Jo J TO 



r 



(2i7rL(ro)Ts) 



Ch{-2j7vs)-l 

2iTTS 



dsdrodr„ 



where the desired signal jhofop ~' r[iVt — n, 1], and [n], 
fc{c) and frifo) are given in ( fTSI ) and Result [T] respectively. 

With single-cell beamforming and no interference coordina- 
tion, the interference at the typical user uo is divided into two 
independent components, the inter-cluster interference /out 
bounded by 7^^^, and the intra-cluster interference denoted 
by /in. The SINR, with perfect CSI beamforming, is given 
by where the desired signal power ||ho|p ^ r[l,A't], and 
/out and /in are independent by the PPP property, since the 
interfering base stations are located in two disjoint areas. 

The probability of coverage with beamforming, with A^t < 
N, is expressed as 



Pc.nic.JVt (Abi Ac, o, T) = 

EN,ro (P [l|h||^ > L(ro)T(/out + hn + 1/SNR) 



(29) 



Conditioned on the value of N, the number of interferers 
inside the cluster is fixed, and the point process corresponding 
to the base stations interferers inside the cluster forms a 
binomial point process (341. The intra-cluster interference term 
/in is thus a result of a binomial point process in the area 
formed by the cluster. For the contribution of /in, we compute 
an upper bound assuming the interference is the aggregation of 
the signal powers of the interfering base stations in the annular 
region with inner radius rp and outer radius tm corresponding 
to the circumscribed circle to the cluster 

The probability of coverage with no interference nulling is 



I ^ /"OO /"CtO /"CtO 

.iVt > < = "] / froiro) fr^irm) fr^ivM) 

[n=Nt ■'° ■'^0 Jro 

/oo 
- OO 



£fa(-2j7r(L(ro)T)-is) - 1 
2jns 



dsdrMdrnidro } , 



where tm is such that its CCDF is bounded by, 1(35 



P[rM > r] > 27rr2e— - (1 + — ^e"'^'' ) 

and the Laplace transform of /^^ is given by 
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Fig. 2. The probability of coverage for increasing SINR tliresliold. Pc.ic 
denotes tlie coverage with per cluster ICIN. Pc.nic denotes the probability 
of coverage with beamforming and no nulling. The average number of base 
stations in the cell is 3 and the pathloss exponent is 4. 



27r((ro + l)(ra.j + l))~°^ 
5^ 5 

^ ((ro + l)2(rM + l)"2Fi(l, 1 - |,2 - |, -(ro + l)-"s) 
-Cm + l)2(ro + l)"2Fi{l, 1 - ^,2 - |, -(rM + 1)-°^)) 
+ (-(^0 + l)(r-M + l)"2Fi(l, 1-^,2- 1 -(ro + 

(r-M + l){ro + l)"2Fi{l, 1 - i,2 - i, -(rM + l)""s)) 

Combining the expressions for Pc,ic,Nt ™d Pcnic.Wt^ we 
obtain the main result on the probability of coverage 
Pc.Nt i^h, Ac, a, T) with interference coordination, for fixed A't 
with thresholding. 

The average throughput for fixed A't follows from the 
probability of coverage analysis similarly to Section HV] 

VII. Simulation Results and Discussion 

We consider a surface comprising on average 100 clusters. 
For random clustering, the density of the cluster base stations 
Ac is varied such that Ac/Ab > 1 to show the benefits of 
increasing cluster sizes. The path-loss exponent is set to a = 4. 

Throughout the numerical evaluation, we compare the per- 
formance of per cluster ICIN with that of unconditional 
beamforming, where each base station beamforms to its own 
users, irrespective of the number of antennas A^t at the base 
stations. With unconditional beamforming, the typical mobile 
user is subject to interference from all the other base stations 
in the network, and the interference power is a shot noise 
process from all the interferers. 

Figure |2] plots the probability of coverage versus the SINR 
threshold T for increasing d^Vt = iVt — iV at the transmitters, 
and for fixed average cluster size Ab/Ac = 3. The probability 
of coverage increases with increasing dAr, , d^^ = 3,d]\[^ ~ 7. 
For low SINR, and for moderate SNR single-cell beam- 
forming outperforms ICIN for d^t = 3. ICIN outperforms 
beamforming as the SINR threshold increases. This can be 
justified by the tradeoff between the decrease in the signal 
power for ICIN with the decrease in the interference power 
at low SINR. At higher SINR values, the system is more 
interference limited, and decreasing the interference when 
using ICIN outperforms boosting the signal power when using 
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SINRnicJVt 



70 1 



^beent,/v('=)(ho)'yo,e\So/e\'^ + T.bgev(-'>{ho)'yo,e\go.e^e 



(l + ro)-"||ho|p 
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SNR 



lout + -fin + 



(28) 
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Fig. 3. The probability of coverage for increasing SINR tliresliold. Pc,ic 
obtained by simulations is compared to the analytical lower bound derived in 
Result 1 . The average number of base stations in the cell is 3 and the pathloss 
exponent is 4. 



beamforming. Moreover, for ICIN, as the clustering is done at 
random, when the base station is at the edge of the cluster, the 
interference from adjacent clusters remains dominant, leading 
to a smaller decrease in the interference power and hence 
lower ICIN coverage gains. For djVt =3, ICIN outperforms 
no nulling in coverage starting at T = 7 dB. This threshold 
decreases with increasing d^Vf For djVt = 7, ICIN outperforms 
no nulling for all the T values of interest. Figure |2] also plots 
the coverage probability performance for fixed lattice grid 
clustering proposed in |[T9l . for d-Nt = 7. We can conclude 
from the comparison that the two models are equivalent in 
terms of coverage performance. 

Figure[3]compares the probability of coverage obtained from 
Monte Carlo simulations with the bound derived in Result 1 
for djVj = 7 and average cluster size equal to 3. Figure[3]shows 
that the lower bound on the probability of coverage exhibits the 
same behavior as the Monte Carlo simulation. It is sufficiently 
accurate, to within 0. 1 in probability, and can provide insights 
on the performance of the clustered coordination system. 

Figure |4] plots the average rate obtained with ICIN as a 
function of average cluster size, for d-Nt =4. It compares the 
average rate with equal bit allocation limited feedback (LF- 
EBA) with and without bias for the desired channel, and the 
average rate with adaptive limited feedback. The average rate 
with no interference nulling is also shown for comparison. 
For a fixed feedback budget at each receiver of Btot ~ 50, 
Figure |4] shows that as the average cluster size increases, the 
average rate with ICIN increases. This is due to the decrease 
in the interference terms as more and more base stations are 
added to the cluster. The same holds for the beamforming 
strategy, as the increase in the number of antennas at each 
base station as a function of the cluster size increases the 
signal term, leading to an increase in SINR. For adaptive 
limited feedback, the average rate increases with the average 
cluster size, similarly to the increase of ICIN with perfect 
CSI. ICIN with adaptive limited feedback outperforms no 
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Fig. 4. The average rate as a function of average cluster size, for A^t = A'^+4. 
r (with no ICIN) also shown for comparison. The feedback budget per mobile 
user is fixed Btot = 50 for all the curves except the LF-EBA curves with 
Btot = {20, 30} indicated in the legend. The performance of beamforming 
(no ICIN) with limited feedback is almost equivalent to the performance with 
perfect CSI. LF-EBA-Bias denotes the equal bit allocation scheme with extra 
bits resulting from rounding to an integer bit value given to the desired user. 
LF-EBA-no Bias denotes the equal bit allocation scheme with the extra bits 
not used. 



ICIN for moderate average cluster sizes up to Ab/Ac = 4. It 
performs almost similarly to no ICIN for larger cluster sizes. 
For limited feedback beamforming with equal bit partitioning 
however, the average rate increases, reaches a maximum at 
Ab/Ac = 2, then decreases. This is justified by the double 
increase in quantization error due to the decrease in the number 
of bits per channel and decrease in the number of bits per 
antenna. Limited feedback with equal bit partitioning and 
biasing towards the desired channel, through allocating the 
extra bits available after rounding the bit allocations to integer 
values, as discussed in Section V, performs better than limited 
feedback with equal bit partitioning and no desired channel 
biasing. The extra bits available at the mobile station from the 
remainder of the division of Btot by N+1 are more judiciously 
used to quantize the desired signal in the equal bit allocation 
with bias. To further illustrate the tradeoff of the average rate 
with equal bit partitioning with the average cluster size, we 
plot the average rate performance of LF-EBA with biasing for 
Btot = {20, 30} in the same figure. For low feedback budget, 
the maximum is achieved at an average cluster size of I. 

Figure |5] plots the mean loss in rate versus the total number 
of bits Btot available at the receiver for Ab/Ac = 3, and 
<i-Nt ~ 5. The figure also shows the analytical upper bound 
derived in Result 3 for comparison. The bound is sufficiently 
tight for all values of Btot- The gap is due to the bounds on the 
quantization errors in the RVQ analysis as well as the Jensen's 
inequality and the Gamma second order moment matching. 
Limited feedback with adaptive bit allocation is shown to 
significantly decrease the mean loss in rate, as compared to 
equal bit allocation. To show the decrease of the mean loss in 
rate for adaptive limited feedback for large values of Btot, we 
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Fig. 5. The mean loss in rate versus the feedback budget per mobile Btot 
is shown. Nt ~ N = 5. and average cluster size is fixed at A^/Ac = 3. 
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Fig. 6. The probability of coverage for increasing SINR threshold. Pc,ic JVt 
obtained by simulations is compared to the analytical lower bound derived in 
Result 1 . The average number of base stations in the cell is 3 and the pathloss 
exponent is 4. 



zoom in on Btot > 30. The mean loss in rate with adaptive 
limited feedback decreases as the number of bits increase, 
albeit at a lower rate than that with equal bit partitioning. 

Figure|6]compares the probability of coverage obtained from 
Monte Carlo Simulations with the bound derived in Section 
VI, for A't = 10 and average cluster size equal to 3. The Figure 
shows that the lower bound on the probability of coverage 
exhibits the same behavior as the Monte Carlo simulation. The 
bound, however is loose for small SINR threshold values. This 
is due to using the circumscribed and inscribed circle to bound 
the interference power for coordination and no coordination. 

Figure [7] plots the average rate obtained with ICIN us- 
ing the thresholding policy in Section VI, as a function 
of average cluster size, for increasing number of antennas 
A^t = {6,8,12}. The achievable rate with unconditional beam- 
forming Tnic.Nt is also shown for comparison. TnicAft does not 
depend on the average cluster size; it is thus almost constant 
for all values of Ab/Ac. Tic.Nt increases with increasing Ab/Ac, 
reaches a maximum, then decreases towards Tnic,Nf This is 
due to the thresholding policy where the feasibility of ICIN 
depends on the number of interferers per cluster When Ab / Ac 
is large compared to the number of antennas, ICIN feasibility 
decreases, and hence single-user beamforming becomes more 
prevalent. The gains from coordination depend on the number 
of antennas at the base stations and the average cluster size. 
For A't ~ 12 for example, as the average cluster size increases. 
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Fig. 7. The average rate r as a function of the average cluster size, for 
increasing A'^t = {6, 8, 12}. r(no ICIN) is also shown for comparison. Tic,jVt 
increases with A^/Ac, reaches a maximum and then decreases. The gains from 
ICIN depend on the number of antennas A'^t at the base stations. 



the average rate increases, reaches a maximum at around 
Ab/Ac =4 and then decreases, until it reaches Taic^N^. As the 
number of antennas increases, ICIN feasibility increases, and 
the average cluster size where ICIN with thresholding attains 
maximum gains increases. The maximum of Tic,Nt depends on 
the average cluster size, the number of antennas A^t and the 
variance of the cluster size. The variance of the cluster size is 
a function of the densities Ab and Ac. As A't increases, Tic,Nt 
increases, and the relative gain from interference coordination 
increases. It reaches 15% for A't = 12 and 8% for A^t = 6. 

VIII. Conclusion 

In this paper, we analyzed intra-cluster interference co- 
ordination for randomly deployed base stations. To cluster 
the randomly deployed base stations, we proposed a random 
clustering strategy that overlays an independent Poisson point 
process on the base stations point process. We assumed that 
the base stations connected to the closest cluster center form 
a cluster. We showed that the average cluster size can be 
optimized with respect to the number of antennas at the base 
stations to provide the maximum gains from ICIN. We further 
analyzed the performance of per cluster ICIN with limited 
feedback CSI. We showed that limited feedback CSI hinders 
the gains of coordination, and results in significant loss in 
rate, when equal bit allocation is used. Adaptive bit alloca- 
tion, optimized as a function of the signal strengths, recover 
the gains of clustered coordination. One takeaway from this 
paper is that analysis of cooperative multi-cell systems, for 
randomly deployed base stations, is possible using hierarchical 
association. Ongoing work includes using the same setup to 
analyze the performance of clustered coordination with CSI 
and user data sharing. It also includes deriving an analytical 
framework for dynamic user-driven clustering, using the same 
framework. 

Appendix A 
Proof of Lemma[T] 

Conditioned on bo being at a distance vq from uo, the 
interference outside a ball of radius — — ri away from 
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uo,the probability of coverage is given by 



Pc,i 



ihSfor > 
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(1 + r-o)' 



'B(uo) ■ 
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{2jTvi{ro)Ts)- 
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\Nt-N 



ra,rm,ri 



1 



-dsdrnidrodri. 



2jTTS 

(a) is derived using the expression in 1361 for the probabiHty 
of coverage using general fading distributions. This expression 
is applicable here since the desired signal power |h|5foP has a 
finite first moment and admits a square integrable density. The 
interference Ir^ is square integrable using the path-loss model 
(l + r)" reminiscent of the pathloss model (l + Ar)" in 1231 . 
The distribution of ri is given by fr{ri) = 2TTXcrie~'^^'^'^'= . 

The aggregate interference considered outside a ball of 
radius r,„ — ro — ri results in an upper bound on the aggregate 
interference due to the fact that less area is excluded from the 
calculations. The inter-cluster interference can also be bounded 
by the interference outside a ball of radius centered at uq, 
corresponding to a slightly larger area to be excluded from the 
calculations than r,,! — ro — ri. This averaged over the spatial 
realizations, results in a good lower bound on the probability 
of coverage as illustrated in Section VII. The lower bound is 
finally given by the expression in Lemma [T] 

Appendix B 
Proof of Remark[T] 

The PMF of N is computed as the Poisson distributed 
number of base stations inside a cluster of size t, such that 
this cluster contains one base station chosen at random from 
the process Ilf,. 

To derive the distribution of the size t of the cluster condi- 
tioned on the availability of one base station chosen at random 
inside the cluster cell, we first need an expression for the 
distribution of the unconditioned size a; of a Voronoi cluster 
As there is no exact result known for the size distribution of 
the Poisson- Voronoi cell, we make use of a two-parameter 
Gamma function fit of the distribution of the normalized cell 
size, derived in 



fxix) 



3.53 



-x'^'^ exp (— 3.5x). 



(30) 



r[3.5] 

Let I be the indicator that a base station chosen at random is 
located inside a cluster cell. The PDF of the size conditioned 
on / = 1 is derived from fx{x) 



fx\i=iit) 



(a) 



fx.i=iit) ^ ¥[I = 1\X = t]f4t) 
F[I = 1] ¥[I = 1] 

3.5*5 



ctfxit) 



r[4.5] 



exp (-3.54), 



(31) 



where c is a constant such that fx\i=i(t)dt = 1, and (a) is 
derived knowing that a randomly chosen base station in lib is 
uniformly distributed in the plane of interest and its probability 
of being inside an area with a given size is proportional to that 



size. The PMF of the number of interferers N inside a cluster 
cell containing one randomly chosen base station is thus 



Pn 



(n) = r¥[N = n\t]fxii=At)dt 
Jo 



(a) 



Jo 



(n)! r[4.5] 



r"'exp(-3.5i)dt 



3.5*-5r(n + 4.5)(Ab/Ac)" 
r(4.5)n!(Ab/Ac + 3.5)"+4-5' 



(32) 



Appendix C 
Proof of Result [3] 

For equal-bit-allocation, the expected residual interference 
is bounded as 



Ne\ (1 + ro,i)-°' 



< E<^ r 12 Ljv+iJiVt-1^ 



(6) 



00 f 
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\Nt-l 



(33) 

where (a) follows by upper bounding the sum of path-loss 
functions from the interfering cells inside each cluster cell 
by the path-loss function from the closest interfering base 
station, ro.i. Because the distances ^ are 1-D PPP with 
intensity TrAb, the random variable TrAbrg i has an exponential 
distribution with parameter 1, ro,i thus has a Rayleigh fading 
distribution ro,i ^ Raylcigh( 1 /-\/27r Ab ) . Therefore 



Ei (l + ro,i)-" 



fro{ro)dxdro. 



(b) follows from expressing the expectation with respect to 
N by its definition as a function of the PMF of N, PN{n). 

To compute E< log2 (cr^ + I, 



.)}. 



, we invoke the Gamma 

approximation of /out- The T{k,9) random variable with the 
same mean and variance as /out with jgg^ffp ~ exp(l) 

(El/ 

has the parameters k and 6 given by fc = varfx) ^^'^ 
9 = f^-^, respectively. The expected value of /out = kO; 
the logarithm of the interference /out — r(fc, 0), with a high 
interference-to-noise ratio (INR) approximation is given by 



(34) 



E| log2(Jo,t)| = log,(2)V.(fc) + \og^{9). 

where tp{k) is the digamma function. 

Appendix D 
Proof of Result |4] 

We consider the terms in the mean loss in rate that cor- 
respond to Bq. Per spatial realization, given TV and ro^e, the 
objective function reduces to 



12 



log2{e)r'^ 



Nt-1, 
"tot - 



(35) 



r ' t , "tot-Up 



'E{/o^t} + l/SN 

We rewrite the objective function as 

Bp 

2 «t-i + 
E{/out} + l/SNR rf "t '1 



= 2 "t-i + C'o2l'=K"t-i) ^ (36) 
and we invoke the arithmetic mean-geometric mean inequality 



2 "t-i + C'o2l'=:i(Nt-i) <2\2 «t-i C'o2 l'=l(Nt-i) 
to find the minimum of the objective function. We solve for 

Bp Bp 

Bo that satisfies the equahty 2 "t~i = Co2 i'ci("t-i) , 

Bo = -(7Vt-l)log2(Co)J|^. 

Replacing Co by its value from (l36l l yields the result for the 
low SNR approximation. 
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